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Summary 

This report reviews aeroelastic analyses for propulsion 
components (propfans, compressors, and turbines) being 
developed and used at the NASA Lewis Research Center. 
These aeroelastic analyses include both structural and aero- 
dynamic models. The structural models include a typical 
section, a beam (with and without disk flexibility), and a 
finite-element blade model (with plate bending elements). 
The aerodynamic models are based on the solution of equa- 
tions ranging from the two-dimensional linear potential 
equation to the three-dimensional Euler equations for 
multibladed configurations. Typical calculated results are 
presented for each aeroelastic model. Suggestions for further 
research are made. Many of the currently available 
aeroelastic models and analysis methods are being 
incorporated in a unified computer program, APPLE 
(Aeroelasticity Program for Propulsion at LEwis). 

Introduction 

An historical perspective of the design and development of 
compressor systems for aircraft gas turbine engines, with 
emphasis on the current research at the NASA Lewis 
Research Center, is presented in reference 1. The author 
noted that further advances in the technology of aircraft 
turbomachinery systems would depend largely on the ability 
to more accurately model and properly account for the 
unsteady flows and their effects in the design process. He 
concluded that solving the problems associated with unsteady 
flow phenomena, such as the loss of stall margin with inlet 
distorted flows, blade flutter within the operating range, and 
premature blade failures due to forced-response excitations, 
have begun to consume more and more of the total develop- 
ment cost of new engine systems. 

Turbine and compressor blade failures due to vibrations 
from unsteady flow can be traced to the infancy of the air- 
craft gas turbine (refs. 2 and 3). The vibrations leading to 
such failures can be stable, as in the case of forced vibrations 
from upstream flow distortions, or they can be unstable, as in 
the case of self-excited vibrations (flutter). In either case, the 
problems are aeroelastic in nature. Addressing these prob- 
lems is important for the development and operation of 
advanced aircraft engines. Data relevant to aeroelastic 
instabilities (flutter) for several types of turbojet and turbofan 


engines are presented in reference 4. These data were 
obtained from a joint NASA/US AF Engine System Research 
program. A synthesis of these data showed that many types 
of flutter can occur near the operating line, triggered by 
different flow conditions (fig. 1(a)). The early reports of 
aeroelastic instability were presented for turbines and 
compressors; however, the newly developed advanced pro- 
pellers (propfans) have also shown a susceptibility to 
aeroelastic instability (refs. 5 and 6). A rotor may also have 
forced-response excitations and fatigue failures from 
dynamic amplification and resonance, which occurs when an 
excitation frequency is near or at a blade natural frequency. A 
Campbell diagram for a rotor blade (fig. 1(b)) can be used to 
identify possible resonant speeds at intersections of exciting 
and natural frequency lines. 

Aeroelasticity is the science that deals with the mutual 
interaction between aerodynamic forces and elastic forces of a 
flexible structure, which in this case is a propulsion component: 
a propeller (propfan), compressor, or turbine. A correct under- 
standing of aeroelastic characteristics is necessary to eliminate 
the problems just described. Aeroelastic analysis methods for 
wings, wing bodies, or complete aircraft (i.e., aeroelastic meth- 
ods related to external flow) have now become fairly well estab- 
lished. Linear methods are available for computations if the flow 
is in either the subsonic or supersonic range. However, for com- 
plex flows containing shock waves, vortices, and flow separa- 
tions, computational methods are needed (refs. 7 and 8). 
Aeroelastic research, in parallel with developments for wing and 
wing bodies, has been in progress to develop aeroelasticity meth- 
ods for propulsion components (i.e., aeroelastic methods related 
to internal flow). A review of this progress is presented in 
references 9 to 11. 

For the last several years, NASA Lewis Research Center has 
been developing aeroelastic analyses for turbomachines and 
propfans. This work has resulted in individual codes with differ- 
ences in the aerodynamic and structural models used. However, 
a single consolidated computer program does not exist for effec- 
tive use by designers and researchers. The availability of numer- 
ous computers from desktop workstations to supercomputers and 
the development of both graphical user interfaces and concurrent 
engineering principles will allow consolidation of all related 
computer codes and aeroelastic analysis methods in a single 
computing system that provides common input and output data 
bases. Such a system, called APPLE (Aeroelasticity Program for 
Propulsion at LEwis), is under development It is the authors’ 
belief that such a system will complement the Numerical 
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(b) Campbell diagram for rotor blade showing possible 
forced-response condition from resonance (denoted 
by circles). 

Figure 1 . — Aeroelastic stability boundaries and forced-response 
conditions for turbomachines and fans. 



Propulsion Systems Simulator (NPSS), also under develop- 
ment at NASA Lewis (ref. 12). 

This paper reviews the different aeroelastic analysis 
methods developed and implemented at NASA Lewis since 
1978 (some of which are documented in ref. 13) for 
turbomachines and propfans. The aerodynamic models vary 
from those based on the linear potential equation to those 
based on nonlinear full potential and Euler equations. For 
aeroelastic analysis, the aerodynamic models were combined 
with structural models that included the typical section, beam 
(with and without disk flexibility), and finite-element blade 
model. Table I summarizes the aerodynamic models (now 
incorporated in APPLE) and their aeroelastic applications. 

Note that the majority of the results presented herein are 
for propfan configurations because they were the primary 
focus of research from 1980 to 1990. However, the analysis 
methods are general and can predict blade flutter and forced- 
response excitations for turbomachines as well. 

This paper is dedicated to the memory of Dr. Krishna Rao 
V. Kaza, who introduced the authors to the field of 
turbomachinery aeroelasticity. 

Multibladed Structures Versus Fixed- 
Wing Aircraft 

For aeroelastic analysis, the blades of a compressor, 
turbine, or propfan (as shown in Fig. 2(a)) are represented as 
a cascade of blades rather than as isolated blades. A 
two-dimensional cascade representation of such a 
multibladed structure is obtained by unwrapping blade sec- 
tions along a stream surface (fig. 2(b)). Phenomena that are 
not encountered with fixed wing surfaces can be attributed to 
the distinctive features of multibladed rotating structures: 
both aerodynamic and structural coupling result from the 
large multiplicity of closely spaced and mutually interfering 
blades; blade attachment and disk flexibility may vary 
greatly; the effect of centrifugal loading is present. Another 
distinctive feature of multibladed structures is that the nature 
of the cascade flow depends on the inlet Mach number and 
the stagger angle; for supersonic flow, the axial component of 
the flow may be subsonic or supersonic, giving rise to differ- 
ent Mach reflection patterns (fig. 3). In addition, two other 
features distinguish the multibladed fluid-structure inter- 
action problems from fixed-wing or isolated airfoil problems: 
(1) structural mistuning and (2) aerodynamic mistuning. 
Structural mistuning refers to slight differences in structural 
properties between the blades of a turbomachine or propfan. 
This mistuning can cause localized mode vibrations in which 
all the energy in the system is concentrated in one or two 
blades, leading to blade failure or loss. Aerodynamic 
mistuning refers to differences in blade-to-blade spacing and 
pitch angles. This mistuning alters the unsteady flow charac- 
teristics in the blade passage. In the present report, only 
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TABLE I.— SUMMARY OF AERODYNAMIC MODELS AND THEIR AEROELASTIC APPLICATIONS 


Aerodynamic model 

Author 

Reference 

Equation 

type 

Aeroelastic application 

Reference 

Model 

type* 

Flow 

I 

Incompressible 

Whitehead, 1960 

30 

2-D 

Kaza and Kielb, 1982 

18 

I 

Subsonic 

Smith, 1973 

31 

2-D 

Kielb and Kaza, 1983 

46 



Rao and Jones, 1975 

32 

2-D 

Kielb and Kaza, 1983 

46 

I 

Transonic 

Surampudi and Adamczyk, 1986 

35 

2-D 

(b) 


I 

Supersonic 

Lane, 1957 

36 

2-D 

Kielb and Ramsey, 1989 

42 


(Supersonic axial) 






I 

Supersonic 

Goldstein, Braun, and Adamczyk, 1977 

33 

2-D 

Busbey, Kaza, and Keith, 1986 

47 


(Subsonic axial) 

Adamczyk and Goldstein, 1978 

34 

2-D 

Adamczyk, Goldstein, and Hartman, 1978 

45 






Kielb and Kaza, 1983 

46 

I 

Subsonic 

Williams. 1985 

43 

3-D 

Kaza, etal., 1987 

50 

n 

Potential 

Verdon and Caspar, 1984 

28 

2-D 

Smith, 1991 

54 

n 

Euler 

Hall and Clark, 1991 

29 

2-D 

(b) 


n 

Navier-Slokes 

(c) 


(c) 

(c) 


m 

Potential 

Kao, 1989 

26 

2-D 

Bakhle, Keith, and Kaza, 1989 

57 



Ku and Williams, 1990 

27 

3-D 

Ku and Williams, 1990 

27 

m 

Euler 

Huff, Swafford, and Reddy, 1991 

24 

2-D 

Reddy, etal., 1991 

62 



Srivastava, 1990 

25 

3-D 

Srivastava, Reddy, and Mehmed, 1992 

65 

m 

Navier-Stokes 

(c) 


(c) 

(c) 



1 Type I, classical linear; Type II, linear, Type III, nonlinear. 
** To be applied. 
c To be developed. 


structural mistuning is considered. Reference 14 provides 
information about systems with aerodynamic mistuning. 

Aeroelastic Formulation and Solution 
Methods 

As mentioned in the Introduction, aeroelasticity involves 
the interaction of structures and aerodynamics. The task of an 
aeroelastic analysis is to combine the formulation of the 
structural dynamic and unsteady aerodynamic models in a 
consistent manner, to solve the resulting aeroelastic model 
for stability and forced vibration, and to interpret those 
results for both qualitative trends and quantitative details 
(fig. 4). A detailed account of the aeroelastic formulation for 
a multibladed structure, or a cascade, is given in reference 15, 
and a brief description follows. 

The formulation for the multibladed structure is begun by 
obtaining the governing equations for A structurally uncoupled 
blades: 

[M]{q} + [*]{*} = {/} (1) 


where [ M ] is the mass matrix, { q ) is the displacement vector, 
[AT] is the stiffness matrix, and {/} is the aerodynamic force 
vector; {/} includes both motion-dependent (self-excited) 


aerodynamic forces which control the stability of the 
aeroelastic system and motion-independent aerodynamic 
forces (due to inlet distortion) which are responsible for the 
forced response. Structural damping can be included in the 
analysis with no additional complication. In the present 
formulation, structural coupling between the blades is 
neglected; hence, the equations for all the blades are coupled 
only through the motion-dependent aerodynamic forces. It 
should be noted further that the Coriolis effects are also 
neglected since most turbomachinery blades are radially 
aligned; studies have shown that the Coriolis effects are also 
negligible for propfan blades. Equations (1) can be solved in 
either the time or frequency domain as shown in the follow- 
ing sections. 

Time Domain Aeroelastic Analysis 

In the time domain aeroelastic analysis, the equations of 
motion (1) are integrated in time starting from the steady 
operating condition and some initial conditions. Thus the 
response, variation of the displacement vector { q } with time, 
is obtained; a growing amplitude indicates the flutter condi- 
tion. The selected steady flow condition determines whether 
the small initial disturbance will grow or decay. The calcula- 
tion is repeated until a flow condition is found for which a 
constant-amplitude motion is obtained that determines the 
flutter boundary. 
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(b) Section along stream surface. 
Figure 2.— A multibladed structure. 


The various methods that can be used to integrate the 
aeroelastic equations in time are given in reference 16. 
Among these, the Newmark method introduces the least 
amount of numerical dissipation and error and thus gives the 
most accurate results. In this method, using constant 
acceleration between two time steps, equations (1) are 
discretized as follows: 

[ M Mr^r + l K M r+Ar {/} T (2a) 

and 


Wt + At= Mr + At Wt 

+ + ( 2b ) 
v H ) 


(a) Subsonic axial velocity. 



(b) Supersonic axial velocity. 

Figure 3. — Cascades in supersonic relative 
inlet flow showing wave reflection patterns. 


(«W {«}, + )({»}, + «U.) <*) 

In equations (2), the subscripts t indicate the time level at 
which the corresponding terms are evaluated. The structural 
dynamic equations are integrated using known aerodynamic 
forces to obtain an updated value of {q}. Then, {q} is used as 
grid motion in the aerodynamic equation solver to obtain new 
aerodynamic forces, which in turn are used to solve for new 
values of {q} and so on. The simultaneous integration of the 
fluid and structure equations is continued until the character- 
istics of the response become clear. If motion-independent 
aerodynamic forces are present, the time domain method 
automatically gives the forced response (variation of the dis- 
placement vector with time); the magnitude of the response 
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Figure 4.— Aeroelastic analysis flowchart. 


where Qa r is the amplitude of the traveling wave of interblade 
phase angle a r , and (0 is the frequency of oscillation. In a 
cascade of N blades, the number of possible interblade phase 
angles is restricted (ref. 17) to 

cr r = 2nr / N r = 0,1,2,-*-,jV- 1 (4) 


Then, the displacement vector for all the blades (assuming 
for clarity that each blade has only one displacement compo- 
nent) can be written as 

{*} = [£]{«„,}*'*” ( 5 ) 

where 


£0,0 

E 0,N-\ 

e n- 1,0 

E N- l.N-l 


can be used to calculate the stress and strain in the blades. 
Because the time domain method does not involve the 
assumption of linearity that is inherent in the frequency 
domain method, it allows the modeling of systems that con- 
tain structural and/or aerodynamic nonlinearities. In addition, 
mistuning effects due to blades with different structural 
properties are accounted for directly and do not need a 
special formulation. 

Frequency Domain Aeroelastic Analysis 

The frequency domain aeroelastic analysis, unlike the time 
domain analysis, is based on the assumption of a linear 
relation between forces and motion. This allows the aerody- 
namic force term {/} in equations (1) to be decomposed into 
motion-dependent j/^Jand motion-independent [f a ] 
contributions as 

{/} = {/<?} + {/*} (3a) 

The frequency domain aeroelastic analysis assumes that it 
is possible to obtain aerodynamic forces for a very specific 
temporal and spatial motion, namely sinusoidal in time and 
fixed interblade phase along the cascade; that is, the n th 
degree of freedom (n th component of the displacement 
vector) of the s ^ blade is written as 

Qns = 2* q o, € ( 3b ) 


F _ ilnsrlN 

Assuming that the blades are undergoing a traveling wave 
motion given by equation (5), the aerodynamic forces J 
and [f a } can be expressed as 

{fa} = ^ 2 « 2 [£][/ ar |^}e' <U( (6a) 


{fa} = ^pb 2 co 2 [E]{f a }e‘ mi (6b) 

where p is the density of the air, b is the semichord, } 
is the complex force coefficient due to blade motion; and 
{/«.) is the complex force coefficient independent of blade 
motion. 

Equations (5) and (6) can be combined with equations (1) 
to yield the governing aeroelastic equations in three ways. 

Traveling wave form . — The governing equations are written 
in terms of interblade phase angle modes by substituting {<7}, 
| /^ J , and [f a } from equations (5) and (6) into equations (1), 
which results in 

-® 2 [£]- 1 [W][£]{^ r } + [Er l [K][E^j 9r } 

= ^%}{K} + npb2c ° 2 { f cr} (7a) 
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Individual blade form . — The governing equations are writ- 
ten in terms of blade displacements by expressing j j in {q s } 
using equation (5) and 

-a) 2 [M]{q} + [*]{?} = npb 2 0) 2 [L]{q} 

+ npb 2 a) 2 {f } (7b) 

where [L] - [£] [la r ] [£*] ' 1 ; {^} is the amplitude of the 
harmonic motion and {/} is the amplitude of the harmonic 
forcing function {/ a }. 

Standing wave mode form. — The governing equations are 
expressed in sine and cosine modes or in structural eigenmodes: 

{?*,} = [£] _1 [^} (8a) 

where 



C s r = cos (2nsr / N) 


S sr = sin (2 nsr I N) 

Substitute equations (6a), (6b), and (8a) into equations (1) 
to obtain the equation for the standing wave form: 

-<o 2 [P)-\m\P^ + [Pr\K]{P)^ 

= npb 2 ( 0 2 [P}- x [L)[P] |jj er j + ap& 2 fi> 2 j^j (8b) 

where 

{/crj = (p r 1 { /a } (&) 

It should be noted that in equations (7a), (7b), and (8b) the 
term e iWt is cancelled from both sides. Also, note that the 
elements of the matrix [E]" 1 are given by 

1 _ 1 „i 2 7Zsr/N 

Ihr r € 

s ' r N 


Once the aerodynamic forces O 0 , ] and [/ 0r }are available 
for all possible values of a r in equation (4) (i.e., all traveling 
wave modes), equation (7) or (8) can be used for stability and 
response prediction. If | f Gr jis set to zero, the equations can 
be cast as an eigenvalue problem and the eigenvalues can be 
used to determine stability. An energy approach can also be 
used to infer aeroelastic stability. If [lo r ] is set to zero, equa- 
tion (7) or (8) results in an algebraic expression from which 
the response can be calculated by simple matrix inversion 
and multiplication. The responses can be summed up to get 
the total response, and they can be used in calculating stress 
or strain on the blades. 


Structural Models 

The various structural models developed thus far are 
shown in figure 5. However, the nature of the formulation 
given in the previous sections allows any other model to be 
included in the aeroelastic analysis provided that the mass 
and stiffness matrices are known. 

Typical Section Model 

The structural model for a cascade consists of a typical 
section with two degrees of freedom (bending and torsion) 
for each blade (see fig. 5(a)). The equations of motion for the 
s th blade are 

m h + S a a + K h h = f h (9a) 

S a h + I a ix + K a (x = f a (9b) 

where m is the airfoil mass, h is the plunging (bending) dis- 
placement, S a is the static unbalance, a is the pitching (tor- 
sion) displacement, 7 a is the moment of inertia, A^and K a are 
the spring constants for plunging and pitching, respectively, 
fh and f a are the aerodynamic loads (including both motion- 
dependent and motion-independent contributions), and the 
dots over the various terms indicate differentiation with 
respect to time. A structural damping term can be included in 
the analysis with very little additional complication. 

The aeroelastic formulation in the frequency domain, 
based on a traveling wave form (eqs. (6)), is detailed in 
reference 18 for both tuned and mistuned cascades; the pro- 
cedures for determining stability and forced response are also 
discussed. These are repeated in the appendix for complete- 
ness. Even though the formulation is presented for a typical- 
section structural model with two degrees of freedom, the 
same formulations can be used with other structural models 
described in the following sections. 
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(d) Finite-element model. 
Figure 5. — Structural models. 
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Beam Model (Ref. 19) 

The beam structural model of each blade consists of a 
straight, slender, twisted, elastic beam with a symmetric 
varying cross section (fig. 5(b)). The elastic, the inertia, and 
the tension axes are each taken to be noncoincident. The 
effect of warping is not explicitly considered. The blade is 
assumed to be rigid in the radial direction. Consequently, the 
radial component of the equations of motion is eliminated. 
The structural model has its basis in the geometric nonlinear 
theory of elasticity in which elongations and shears are 
negligible compared with unity, and the squares of the 
derivatives of the deformation of the elastic axis are negli- 
gible compared with the squares of the bending slopes. This 
level of the geometric nonlinear theory of elasticity is 
required to derive a set of linear coupled bending-torsion 
equations of motion. 

The blade degrees of freedom consist of bending in the 
plane of rotation u , bending in the plane perpendicular to the 
plane of rotation w , and pitching (torsion) about the elastic 
axis, a, as shown in figure 5(b). In deriving the equations of 
motion, three coordinate systems are used (fig. 5(b)). The 
Xq, is the hub fixed system and rotates about the Zq 

axis with an angular velocity Q; x, y, z is the blade fixed 
system at an arbitrary point 0 on the elastic axis, and X 3 , y 3 , 
Z 3 is the blade fixed system at O' after deformation, obtained 
by rotating xyz by a. The equations of motion are derived 
using the extended Hamilton’s principle. By substituting 
expressions for strain energy, kinetic energy, and virtual work 
of aerodynamic forces and by performing the required opera- 
tions, the governing equations are obtained. The primary 
motions are expanded in terms of the generalized coordinates 
{ q }, associated with the nonrotating uncoupled beam modes 
in pure bending and torsion, to obtain the governing equa- 
tions in the form shown in equations (1). The elements of [A/] 
and [A"] are defined in reference 18. Thus, formulations for 
both tuned and mistuned systems can be written in a form 
similar to those presented for the typical-section model. 

Bladed-Disk Model (Ref. 20) 

The bladed-disk system is idealized as a constant-thick- 
ness, uniform circular disk to which the blades are attached; 
only the disk out-of-plane bending motion uj is considered 
(fig. 5(c)). With the exception that the warping of the blade 
cross sections is explicitly considered, the blade model is the 
same as that discussed in reference 19 and explained in the 
section Beam Model. The equations are written in traveling 
wave form for the blades (eqs. (7)) and in standing wave 
form for the disk (eqs. ( 8 )). Continuity and slope continuity 
are enforced at the blade-disk junction. Hamilton’s principle 
is again used to obtain the matrices [Af] and [A"]. 


Finite-Element Model (Ref. 21) 

Although the models just described help in understanding 
the basic aeroelastic problem, a finite-element model 
(fig. 5(d)), is often required to accurately represent more 
complex blade structures such as propfan blades. Because 
propfan blades are thin and flexible, deflections due to 
centrifugal and aerodynamic loads are large. Hence, the 
aeroelastic problem is inherently nonlinear, requiring the 
geometric nonlinear theory of elasticity. The blades of a 
propfan have a low aspect ratio and they behave as platelike 
structures having chord flexibility and requiring a 
three-dimensional structural model. They have large sweep 
and twist, which couples blade bending and torsional 
motions. These factors require a finite-element structural 
model that also accounts for centrifugal softening and 
stiffening effects. Using finite-element models, the governing 
equations can be written again as shown in equations (1). The 
matrices [Af] and [A'] are obtained from the finite-element 
analysis. Again, these can be used in the aeroelastic 
formulation with appropriate aerodynamic models. For most 
of the studies at Lewis, a NASTRAN finite-element analysis 
is used. A review of the available finite elements in COSMIC/ 
NASTRAN and MSC/NASTRAN, along with their 
capabilities and limitations for modeling flexible rotating 
blades, is given in reference 22; reference 23 describes the 
geometric nonlinear analysis using MSC/NASTRAN for 
frequency and mode shape calculation of blades under 
centrifugal loading. 

Aerodynamic Models 

The various aerodynamic models used at Lewis are 
presented in table I. Nonlinear models (refs. 24 to 27), 
referred to as Type III models, include the effects of the blade 
shape and thickness. In addition, the nonlinear models are not 
restricted to small-amplitude blade oscillations so the 
assumption of sinusoidal motion is not required. These 
models are valid for arbitrary motion but require more 
computing time. Two types of linear unsteady cascade 
analyses have been developed on the assumption that the 
unsteady disturbance is small. One analysis (Type II) linear- 
izes about a nonuniform, deflected mean flow (refs. 28 
and 29); the other (Type I) linearizes about a uniform, 
undeflected mean flow (refs. 30 to 36). In the present report, 
Type II models are referred to as linear models and Type I 
models as classical linear models. Since these linear models 
are based on a small perturbation analysis, nonlinearities 
related to the amplitude of motion are not modeled in 
either analysis. 
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Type III (Nonlinear) Models 

The nonlinear models require the numerical solution of the 
equations of fluid dynamics. In reality, both steady and 
unsteady flows in turbomachines and propfans are extremely 
complicated. The fluid is viscous and heat conducting and is 
most accurately described by the Navier-Stokes equations. 
However, if the Reynolds number is sufficiently high and 
separation does not occur, then the viscous and heat-transfer 
effects are confined to narrow regions near the airfoil sur- 
faces and the wakes. Under these circumstances, the Euler 
equations provide a good approximation of the behavior of 
the flow. The three-dimensional Euler equations in a 
Cartesian frame are given as 


dQ dE dF dG . 

dt dx dy dz 

(10) 

where 


Q = [p,pu,puv,pw,e] T 

(11a) 

E = [pu,pu 2 + p,puv,puw,u(e + p)j 

(lib) 

F = [pv, puv, pv 2 + p, pvw, v(e + p)] 

(lie) 

T* 

G = [pw, puw, pvw, pw 2 + p, w(e + p)j 

(lid) 

and p is the fluid density, u, v, and w are the Cartesian veloc- 
ity components, e is the energy per unit volume, and p is the 
fluid pressure. Appropriate boundary conditions can be 
specified to complete this formulation. 

Generally, equations (10) in conservative differential form 
are transformed from a Cartesian reference frame to a 


time-dependent, body-fitted curvilinear reference frame for 
numerical solution. The transformed two-dimensional Euler 
equations are solved using the flux-difference-splitting (FDS) 
scheme in reference 24. The three-dimensional Euler 
equations are solved using the alternating-direction-implicit 
(ADI) scheme for propfan configurations in reference 25. 

For irrotational flow, the full potential equation is obtained 
from the Euler equations: 

dp [ djpu) [ djpv) | d{pw) _ Q ( 12 j 

dt dx dy dz 


where 



In the preceding equations, <p is the velocity potential; a«„ 
and Mac are the density, sonic velocity, and Mach number, 
respectively, at the reference condition. Kao (ref. 26) solved 
the two-dimensional version of the full potential equation in 
finite-volume form by using the Newton iteration method for 
a cascade of blades. The three-dimensional equations for 
propfan configurations have been solved by Ku and 
Williams (ref. 27). 

The nonlinear unsteady aerodynamic models allow either a 
frequency domain or a time domain aeroelastic analysis to be 
performed using the same aeroelastic solver. Since the 
frequency domain flutter analysis is linear in nature, the use 
of a nonlinear unsteady aerodynamic model in a frequency 
domain flutter calculation may seem inconsistent, but it is 
justified in the absence of a corresponding linear aerody- 
namic solver. The unsteady harmonic force coefficients 
required in the frequency domain analysis can be obtained by 
harmonically oscillating the blades and by Fourier- 
decomposing the resulting forces. However, this process is 
time consuming in comparison with those of the Type I and 
Type II models in which a harmonic variation of the flow 
variables is assumed, thus removing the time dependency. In 
order to make the time-accurate nonlinear codes more 
computationally efficient, two methods have been developed 
and implemented in the codes: the Influence Coefficient and 
Pulse Response methods (ref. 37). These methods are valid 
for small-amplitude blade oscillations for which the unsteady 
flow field is linearly dependent on the amplitude. Other 
methods of reducing computational time such as parallel 
processing (ref. 38) and reduced-order modeling for flutter 
analysis (ref. 39) are also pursued at Lewis. These methods 
reduce the. computational time required and allow the nonlin- 
ear codes to be used efficiently for a linear frequency domain 
analysis while the flexibility of a time domain analysis is 
maintained. 

In addition, viscous terms can be easily included in the 
nonlinear Euler equations and can be used for accurate 
prediction of stall and choke flutter (see fig. 1(a)), which 
involve nonlinear behavior. Earlier methods have used 
empirical models to investigate these types of flutter 
(refs. 40 and 41). 
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Type II (Linear) Models 

Type II models use full nonlinear equations for the steady 
solution but use linearized equations for the unsteady 
solution (refs. 28 and 29). Thus, the effect of airfoil shape 
and angle of attack can be included in the unsteady aerody- 
namic and aeroelastic calculations. 

The model developed by Verdon and Caspar (ref. 28) is 
based on the two-dimensional potential equation. The poten- 
tial function is expanded in a perturbation series as 

<p(x,y,t) = <fc)(x,y) + £<t>\{x,y,t) + ... (13) 

where € is a small parameter. Substituting back into the 
governing equation (eq. (12)) and neglecting terms of 0(e 2 ), a 
nonlinear equation in <po for steady flow and a linear variable- 
coefficient equation in <f>\ for unsteady flow are obtained. The 
equation for steady flow is the steady nonlinear (full) 
potential equation; this equation is solved in order to obtain 
the steady flow field, which in turn is used to calculate the 
coefficients in the unsteady equation. A further simplification 
is made by assuming that the flow field varies harmonically 
with time. This assumption allows a single interblade passage 
of the cascade to be used for calculations with any interblade 
phase angle and also removes the time dependence from the 
formulation. Hall and Clark (ref. 29) developed a similar 
model based on two-dimensional Euler equations. 
Three-dimensional models are in a developmental stage. 

Type I (Classical Linear) Models 

In Type I models, the unsteady potential flow is linearized 
about a uniform mean flow; that is, 0 o is linear. Thus, the 
effects of blade shape (airfoil thickness and camber) are 
neglected entirely (refs. 30 to 36). These models can be 
applied only to cases in which the effects of blade geometry 
can be neglected. Whitehead, Smith, and Rao and Jones 
(refs. 30 to 32) solved the two-dimensional linear potential 


equation for subsonic flow by using different numerical tech- 
niques. Goldstein, Braun, and Adamczyk (refs. 33 and 34) 
solved the linear potential equation for supersonic flow with 
a subsonic leading edge locus. Surampudi and Adamczyk 
(ref. 35) solved the equations for nominally sonic flows. The 
formulation by Lane (ref. 36) is coded into a program by 
Kielb and Ramsey (ref. 42) for solving supersonic flow with 
a supersonic axial component. A formulation based on linear 
compressible small-disturbance theory was developed for a 
three-dimensional rotating multibladed geometry (ref. 43). In 
this work, an integral equation relating the normal velocity to 
the load distribution on the blade is derived. The integral 
equation is solved by dividing the blade into a finite number 
of elements (panels) on each of which the load is constant. 
The loads on all the elements are determined simultaneously 
by requiring that the normal induced velocity be equal to the 
specified value at the control point. This formulation was 
recently extended to single-rotation propfans with 
ducts (ref. 44). 

Results and Discussion 

The structural and aerodynamic models described in the 
previous sections were applied to investigate the aeroelastic 
stability and response of wind tunnel propfan models. The 
following selected results are a comparison of analyses and 
experimental measurements. It should be noted that the 
analyses are limited to isolated blade rows or to single- 
rotation (SR) propfans. Table II presents the propfan configu- 
rations that were referenced in this report. 

Type I (Classical Linear) Aerodynamic Models 

The earliest aeroelastic model was the typical-section 
structural dynamic modeI>4thj»lunging and pitching 
motions of each blade along with a tvto-dimensional linear 
cascade aerodynamic model. This model was used as a 
research tool to understand the physics of cascade effects. 


TABLE C.— PROPFAN CONFIGURATIONS 


Model 

Configuration 

Tip 

geometric 

sweep, 

deg 

Number 
of blades 
in full 
rotor, 

N 

Material 

Reference 
diameter 
m (in.) 

Ratio of 
hub to 
blade 
diameter 

SR3 

Single rotation 

45 

8 

Titanium 

0.622 (24.5) 

0.239 

SR3CX2 



45 

8 

Composite 

.622 (24.5) 

.239 

SR3C-3 



45 

8 

Composite 

.622 (24.5) 

.239 

SR5 



60 

10 

Titanium 

.622 (24.5) 

.235 

SR7L 



41 

8 

Composite shell/ 

2.74(108) 

.232 






titanium spar 



SR7A 



41 

8 





F21 

Counter rotation 

45 

13 



.622 (24.5) 

.232 


(forward rotor) 





.617 (24.3) 

.430 
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Researchers who investigated flutter in tuned rotors in 
supersonic flow with a subsonic leading edge locus used 
an energy approach to infer the stability (ref. 45). The flutter 
analysis was applied to compressors at both high and 
low backpressure conditions by using the theories of 
references 33 and 34, respectively. The calculated stability 
trends were consistent with experimental rig observations. In 
reference 18, the typical-section model was also used to 
understand the effects of mistuning; an eigenvalue approach 
was employed to determine stability. The results obtained 
using this aeroelastic model showed that bending-torsion 
coupling has a significant effect on cascade flutter. Also 
noted was that frequency mistuning had a beneficial effect 
on suppressing flutter in all the flow regimes considered 
(refs. 18, 42, 46, and 47) and that it had either a beneficial or 
an adverse effect on forced response. These results were later 
used as benchmarks for checking more complicated models. 

Subsequently, the beam structural model was used along 
with a two-dimensional aerodynamic model applied in a 
stripwise fashion (ref. 19). Disk flexibility was included in 
the aeroelastic formulation in reference 20. Parametric 
studies indicated that the effect of frequency mistuning on 
flutter was still beneficial even in the presence of structural 
coupling between blades due to disk flexibility. Also, it was 
found that blade pretwist introduces strong coupling between 
the disk bending and blade chordwise motions. 

The beam model was used to predict and correlate the flutter 
speed of the SR5 propfan wind tunnel model and to clarify the 



-# Experimental flutter 



Figure 7. — Comparison of beam model flutter 
boundary and experimental flutter points for 
SR5 propfan (ref. 5). 


mechanism of the flutter phenomenon (ref. 5). This beam model 
was modified to account for blade sweep in an approximate 
manner as shown in figure 6. The two-dimensional aerodynamic 
models were also modified to account for sweep by using 
similarity laws. The disk was assumed to be rigid. Figure 7 
shows the calculated flutter boundary in comparison with 
experimental data. The measured and calculated flutter 
boundary trends are seen to be in good agreement, as were 
the flutter frequencies. Also, the flutter interblade phase 
angles agreed well with the calculations. However, because 
the analytical results required judgement in selecting an 
effective blade sweep and a blade elastic axis position, a 
three-dimensional finite-element structural model was devel- 
oped to circumvent making such a decision. 

A finite-element structural model and a three-dimensional 
aerodynamic model were applied to correlate the flutter 
boundary of the SR3CX2 propfan (ref. 21). A NASTRAN 
finite-element model was used to obtain the mode shapes and 
frequencies of the propfan. The three-dimensional linear 
aerodynamic model of reference 43 was used. The calculated 
shapes and frequencies of the first two natural modes are 
shown in figure 8(a). The mode shapes include the effect of 
centrifugal load. The following findings were made: the first 
mode was primarily bending and the second mode primarily 
torsion; centrifugal loads increased blade twist and aerody- 
namic loads decreased twist; the combined loads resulted in a 
net increase in twist, and this twist increased with rotational 
speed. In comparison with measured flutter data, the calcu- 
lated flutter results are presented in figure 8(b) for four and 
eight blades. The difference between the calculated and 
measured flutter Mach numbers was greater for the 
four-blade case than for the eight-blade case, which implies 
that the theory may not be accurate for large blade spacing. 
Calculated and measured interblade phase angles also 
compared well; however, the calculated flutter frequencies 


ll 



0 




Mode 1 2 

Frequency, Hz 200 399 

(a) Calculated mode shapes and frequencies. 



Free-stream Mach number, 

(b) Measured and calculated flutter boundaries. 


Figure 8. — SR3CX2 propfan mode shapes, frequencies, and flutter boundaries at blade setting 
angle p ■ 61 .6° (ref. 21 ). 


were about 8 percent higher than the measured. Calculations 
were also performed using the three-dimensional structural 
model with two-dimensional aerodynamic models applied in 
a stripwise (quasi-three-dimensional) manner for the 
eight-blade configuration. The correlation varied from 
poor to good as shown in figure 9. Recently, this 
quasi-three-dimensional aeroelastic analysis program was 
used to design flutter-free cruise missile propfans for a joint 
Navy/Air Force/NASA project (ref. 48). A similar 
quasi-three-dimensional aeroelastic analysis was also 
implemented in a NASTRAN analysis (ref. 49). 

A three-dimensional aeroelastic model was used to 
perform the calculations for a propfan having alternate blade 
mistuning (alternate blades have different structural proper- 
ties). These calculations are presented in reference 50. The 
propfan had four SR3CX2 blades and four SR3C-3 blades 
which differ in natural frequencies and mode shapes. The 
measured and calculated flutter boundaries are shown 
in figure 10. The overall agreement between theory and 
experiment is good. 


The three-dimensional unsteady aerodynamic model 
(ref. 43) has the capability of performing a modal 
forced-response vibration analysis including structural and 
aerodynamic mistuning; the excitation is aerodynamic in 
origin. August and Kaza (ref. 51) analyzed an aeroelastically 
scaled model, SR7A, which was used to simulate a large- 
scale propfan, SR7L; the performance, vibration, forced 
response and flutter of both propfans were studied. They con- 
cluded that the aeroelastic model accurately simulated the 
prototype fan. The measured and calculated one-per- 
revolution, forced vibratory stress amplitudes for the root 
region of the SR5 propfan blade are shown in figure 1 1 (from 
ref. 52); the calculations are based on the three-dimensional 
unsteady aerodynamic model. Good correlation between 
experiment and analysis is seen; the correlation for the tip 
region (not shown here) is not as good. The flutter predic- 
tions for SR5 from reference 13 are shown in figure 12. The 
maximum difference between theory and experiment was 
16.7 percent at 6800 rpm, and the numerical predictions are 
seen to be unconservative. It should be noted that, even 
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Figure 9. — Flutter results from quasi-three dimensional and 
three-dimensional aeroelastic models for SR3CX2 propfan 
(ref. 21). 


Rotor 

O SR3CX2 
□ SR3CX2/SR3C-3 

Solid symbols denote 
measured data 


Open symbols denote 
values calculated with 
measured second-mode 



Figure 10.— Flutter boundaries for mistuned propfan (ref. 50). 


though both stress and flutter predictions fall within 
20 percent of the experimental data, the stress prediction will 
result in a safe design whereas the flutter prediction will 
result in an unsafe design. The observed differences indicate 



Inflow angle, deg 

Figure 1 1 .—Vibratory stress amplitudes 
for SR5 propfan based on three- 
dimensional aerodynamic model (ref. 
52); Mach number, 0.8; rotational 
speed, 5500 rpm. 


0.85 0.75 0.71 Tunnel 

Mach 

number 



Figure 1 2. — Flutter boundaries for SR5 propfan (ref. 1 3) at blade 
setting angle Po.75R = 69.3°. 


the need for more accurate modeling in the tip region and for 
the modeling of transonic flow, which would require 
advanced aerodynamic models (Types II and III). 
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Imaginary part of eigenvalue, ?/w, 



Mode 

Frequency, 

Hz 

Bending 

4 748 

Edgewise 

9 950 

Torsion 

16 580 


(a) Finite-element model and natural frequencies; 10014 nodes; 
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(b) Root locus plot from two-dimensional linearized potential 
solver for vibration in second (edgewise) mode, zero struc- 
tural damping. 


Figure 13.— Space shuttle main engine (SSME) high- pressure- 
oxygen turbopump (HPOTP) turbine aeroeiastic analysis 
(ref. 54). 


Type II (Linear) Aerodynamic Models 

As can be seen from table I, the only aeroeiastic applica- 
tion of the Type II aerodynamic model is the one based on the 
two-dimensional linearized potential solver (ref. 28). This 
model improves upon the classical linear aerodynamic mod- 
els used in the previous section in the following two ways. 


The effects of airfoil shape and angle of attack are included 
in the aerodynamic model, which has the capability of treat- 
ing arbitrary modes including those with chordwise flexibil- 
ity. A typical-section structural model was used to study the 
effect of steady aerodynamic loading on the flutter of a com- 
pressor cascade (ref. 53). The study showed that the neglect 
of steady aerodynamic loading in flutter calculations could 
result in nonconservative estimates of the flutter boundary. 
This aerodynamic model was also applied to investigate the 
stability and forced response of the space shuttle main engine 
(SSME) high-pressure-oxygen turbopump (HPOTP) turbine 
(ref. 54). The blades in the first stage of the turbine experi- 
enced frequent cracking in the shank region. The aeroeiastic 
analysis was formulated in modal form, and the aerodynamic 
model was applied in strip wise fashion. The MSC/NASTRAN 
finite-element program was used to perform the structural 
dynamic analysis of the HPOTP turbine blade. The model 
consisted of 10 014 nodes and 7758 solid brick elements, as 
shown in figure 13(a) along with calculated natural frequen- 
cies. The span of the blade was divided into eight strips at 
each of which a two-dimensional linearized aerodynamic 
analysis was then applied. The flutter analysis indicated that 
the HPOTP turbine blades experienced very low aerody- 
namic damping in the first four vibrational modes. The sec- 
ond mode (first edgewise mode) was unstable (fig. 13(b)). An 
addition of 1 percent damping made the system stable. A 
mistuning analysis of the HPOTP blade (ref. 55) showed that, 
when small mistuning (of the order present in actual rotors) 
was introduced, the aeroeiastic modes become localized to a 
few blades, possibly leading to blade failure. Some prelimi- 
nary calculations of the HPOTP blade forced response due to 
gust and cooling jet were presented in reference 56. 

Type III (Nonlinear) Aerodynamic Models 

The nonlinear models were developed to investigate com- 
plicated flows involving shocks and nonlinearities in the 
flow. However, these models were validated in various limit- 
ing cases by comparison with linear and classical linear 
theory. The two-dimensional full potential equation was used 
to predict and validate flutter calculations for selected cases 
in references 57 to 60. In reference 58, a typical-section 
structural model was used to simulate the SR5 propfan. The 
aeroeiastic equations were solved both in time and frequency 
domains. In the calculations, the cascade had five blades and 
the Mach number at the inlet was 0.7. Figure 14(a) shows the 
root locus plot with all 10 eigenvalues, 2 for each of the 
5 blades in the cascade, at a reduced frequency k c = 0.225. 
The eigenvalues fall into two groups. In the group at the 
lower frequency, the most unstable phase angle is 288°. 

A time domain solution was also obtained by the integra- 
tion of the equations of motion for each blade. Figure 14(b) 
shows the variation of the center blade pitching displacement 
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(a) Frequency domain solution for k c * 0.225 (k c = 



(b) Time domain solution for two values of reduced 
velocity V* = M 00 a 00 /bw a . 

Figure 14.— Flutter calculations for simulated SR5 propfan 
typical-section model from two-dimensional fuH potential 
solver; Mach number, 0.7 (ref. 58). 

with time for stable and unstable conditions. The flutter 
velocity, frequency, and interblade phase angle showed good 
agreement with the corresponding results from the frequency 
domain method, thus verifying both the time and frequency 
domain aeroelastic analysis methods. The variation of the 
flutter velocity with mach number to predict the transonic dip 
is presented in reference 59. It was found that the transonic 
dip in cascades is similar to that for isolated airfoils (ref. 61), 
except that it occurs at significantly lower Mach numbers 
(see fig. 15). In reference 27, the three-dimensional full 
potential solver was used to analyze SR3CX2 propfan, for 
in-phase motions (a = 0°) both in the time and frequency 
domains. The calculations and the experiments showed the 
blades to be stable. 

The two-dimensional Euler solver was applied to calculate 
the flutter behavior of selected typical-section models in 
references 62 and 63. A calculation of the flutter condition 
for the SR5 propfan model produced the same results as 
those calculated by the full potential solver; these results are 



Figure 1 5. — Transonic dip for cascade with NACA 64A01 0 
airfoil using full potential aeroelastic solver (denoted as 
''present") with comparisons to classical linear theory and 
isolated airfoil theory (ref. 59). 

summarized in table III. The ability to predict transonic flow 
is demonstrated in figure 16 for an example NACA 0006 cas- 
cade (ref. 28). This figure shows the steady pressure distribu- 
tion and a comparison with a full potential equation solution. 
The Euler equations predict the shock location slightly down- 
stream from that predicted by the full potential solver. The 
unsteady pressure distribution is compared with a classical 
linear analysis in figure 17. This comparison shows the effect 
of airfoil shape on the unsteady loading and clearly indicates 
that classical linear theories which neglect the effects of 
airfoil shape are inadequate for transonic flow calculations. 

The three-dimensional Euler equations were used in 
reference 64 to study the effect of structural flexibility on the 
performance of SR7L propfan, by coupling the Euler aerody- 
namic solver with a NASTRAN finite-element analysis. It 
was concluded that the customary way of adjusting the blade 
setting angle by rigid-body rotation does not correctly simu- 
late the actual blade shape during operation. The 


TABLE III. — COMPARISON OF FLUTTER BOUND- 
ARIES CALCULATED FROM FULL POTENTIAL 
AND EULER SOLVERS* 


Parameter 

Full 

potential 

Euler 

Interblade phase angle, (T, deg 
Reduced frequency, k ( - ox fM jct^ 
Frequency ratio, 0) /(t) a 
Reduced velocity, Vj = M / ba> n 

288 

0.225 

.61 

5.43 

288 

0.222 

.67 

6.10 


a References 58 and <S3. 
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Distance along chord 

Figure 1 6. — Steady Mach number distribution for transonic 
cascade; gap-to-chord ratio, 1 .0; stagger angle, 45°; relative 
Mach number, 0.8; incidence angle, 13° (ref. 63); full potential 
solution is denoted as Verdon (28) in the key. 



Figure 1 7. — Unsteady pressure distribution for transonic cas- 
cade; pitching about midchord; parameters as In figure 16 
(a); reduced frequency (based on semichord) of 0.5, interblade 
phase angle of 0°, amplitude of oscillation of 2° (ref. 63); 
linear theory results are denoted as Smith (31 ) in the key. 

three-dimensional Euler equations were also solved in refer- 
ence 65 to verify the static stability characteristics of two 
forward-swept propfan blades. The analysis confirmed the 
static stability behavior of the blades that were designed with 
uncoupled aerostructural analysis. The study also confirmed 


that a forward-swept blade can be designed to be statically 
stable by tailoring the structural properties of the blades using 
composites. 

Flutter calculations have been performed for the F21 
propfan using both the three-dimensional full potential and 
Euler solvers. Neither of the codes predicted the observed 
flutter of the F21 propfan in the in-phase (zero interblade 
phase angle) mode. The presence of a leading edge vortex at 
the experimentally observed flutter condition is considered 
significant in this case. The full potential formulation does 
not allow the modeling of such a vortex. Preliminary Euler 
calculations show the presence of a leading edge vortex. 
Further verification using finer grid spacings is required 
before accurate flutter predictions can be made including the 
influence of the vortex. Both the full potential and Euler 
aeroelastic solvers, which were previously restricted to 
in-phase blade motions, have been extended to allow the 
modeling of multiple blade passages with independent blade 
motions. This will allow further validation of these solvers 
using experimental data not limited to the in-phase 
flutter mode. 

APPLE 

The structural and aerodynamic models and the time and 
frequency domain solution methods that have been 
developed and validated are being incorporated in a 
user-oriented program with a common data base for input and 
output. The acronym for this program is APPLE 
(Aeroelasticity Program for Propulsion at LEwis). A 
flowchart of this system is shown in figure 18. A graphical 
user interface will provide the researcher or analyst with an 
easy way of selecting the desired aeroelastic model suitable 



Figure 18. — Various options available in APPLE system. 
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for his/her purpose. This program, when completed, will run 
on a number of distributed workstations using parallel 
processing, a graphical user interface, and concurrent 
engineering methods. The APPLE system, when developed, 
can be directly integrated in the Numerical Propulsion 
Systems Simulator (NPSS) project being developed at 
NASA Lewis. 

Summary of Results 

A review of the aeroelastic models available for propulsion 
components at the NASA Lewis Research Center was pre- 
sented. It is evident from this material that a number of 
aeroelastic models and methods exist for investigating flutter 
and the forced response of isolated blade rows and 
single-rotation propfans. The models developed thus far were 
valid, though restricted to single-rotation propfans, for a 
variety of operating conditions. 

Recommendations 

The following suggestions are made for future research: 

(1) A comprehensive study is required to verify and 
evaluate the transonic modeling capability of nonlinear 
aeroelastic analyses for transonic flutter prediction. 

(2) The development and aeroelastic application of 
three-dimensional linearized (Type II) aerodynamic models, 
based on both the potential and Euler equations, is required. 

(3) All the models surveyed in this paper were restricted to in- 
viscid flows in which viscous effects were neglected. For appli- 
cations in which viscous effects are expected to be significant, 
viscous flow models need to be developed. 

(4) Methods to predict aerodynamic forces due to inlet distor- 
tions, incident pressure waves, and other forcing functions need 
to be developed. 

(5) Stall flutter and choke flutter remain significant problems 
in the design and operation of axial compressors and fans. In the 
past, empirical models have been used to analyze these types of 
flutter. Type III viscous and transonic flow models need to be 
developed to predict these types of flutter without using empiri- 
cism. 

(6) Results presented in this report revealed that flutter usually 
occurred in the first two or three modes (i.e., low-frequency 
modes). The structural analysis methods presented earlier 
showed good correlation for frequencies and mode shapes in 
these low-frequency modes. However, during experiments with 
low-aspect-ratio compressor blades made of composite materi- 
als, it was found that flutter can occur at higher frequency modes. 


In these modes, the largest vibratory motion is usually at the tip. 
The high frequency vibration leads to high cycle fatigue and loss 
of the tip of the blade. Analysis methods that predict accurate fre- 
quencies and mode shapes at higher modes should be developed. 

(7) In order to meet aeroelastic and forced-response 
constraints, blades of large jet engine fan rotors are often 
stiffened by one or more sets of part-span shrouds. A major unre- 
solved problem in the dynamic analysis of shrouded rotors is the 
determination of the appropriate boundary conditions at the 
blade-shroud interfaces. Nonlinearities are introduced in both the 
structural and the damping operators because of the nature of dry 
friction damping and the dependence of boundary conditions on 
slip and blade displacements. These nonlinearities can give rise 
to subcritical instabilities and jump phenomena which have a 
bearing on flutter predictions. Earlier analysis methods used a 
linearized approach to treat these nonlinearities. With the devel- 
opment of time domain aeroelastic methods, it is now possible to 
treat nonlinearities more accurately; this approach should be 
pursued. 

(8) Mechanical damping in turbomachines occurs during rub- 
bing at mating interfaces: shroud, root, blade-to-blade or blade- 
to-disk dampers. Modeling or characterizing material damping 
has also been a problem for a long time. More accurate modeling 
of mechanical and material damping is needed for making better 
aeroelastic predictions. 

(9) System analyses are required for making reliable blade 
and aeroelastic response predictions. The analyses must include 
multicomponent coupling, such as blade-disk and disk-rotor in- 
teractions and support system flexibilities. Substructuring algo- 
rithms can provide the necessary aeroelastic information while 
maintaining accurate overall system responses. 

(10) Propulsion research at NASA Lewis is now focusing on 
advanced ducted fans and turbines. The ultra-high-bypass (UHB) 
engine concept with ducted, wide-chord fan blades is an example 
of such trends in engine design. Existing analyses methods need 
to be modified for these ducted turbo machinery configurations. 

(11) Advanced engines require aeroelastic analysis for mul- 
tiple blade rows, such as rotor-stator stages and counter-rotating 
propfans. Aeroelastic methods that account for unsteady aerody- 
namic interactions and structural coupling for these configura- 
tions should be developed. 

(12) The higher flight Mach numbers planned for the High- 
Speed Civil Transport (HSCT) will require that thermal analysis 
be coupled with aeroelastic analysis. Also, acoustic should be 
developed. Aeroelastic tailoring methods with flutter and blade 
strain constraints should also be developed for efficient light- 
weight propulsion components. 

(13) Reduction in computational time is a requirement if the 
models developed are to be useful for design and optimization 
applications. Methods that use parallel processing and reduced- 
order models should be vigorously pursued. 
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Appendix-Formulation and Procedure for Aeroelastic Stability Analysis of Cascades 


For completeness, the aeroelastic formulation in the fre- 
quency domain (from ref. 18) is presented for the stability 
analysis of tuned and mistuned cascades: The procedure for 
predicing the flutter boundary from these equations is also 
presented. The typical-section structural model is used as an 
example in presenting the equations. The extension to more 
refined structural and aerodynamic models is straightforward. 

Aeroelastic Analysis for Mistuned Cascades 

The following analysis is presented for a generally 
mistuned cascade in which each blade may have different 
structural properties. The analysis for the special case of a 
tuned cascade in which all blades are identical is presented in 
the subsequent section. The approach followed assumes that 
the structure is vibrating in an aeroelastic mode (interblade 
phase angle mode) with a motion that is a harmonic function 
of time. The frequency of oscillation is permitted to take on 
complex values, thus allowing decaying-, growing-, or 
constant-amplitude oscillations. The aerodynamic forces cor- 
responding to constant-amplitude harmonic oscillations are 
inserted in the equations of motion to formulate a complex 
eigenvalue problem. The eigenvalues are generally complex 
quantities, and therefore a complex frequency is obtained. 
The imaginary part of the complex frequency represents the 
damping ratio and thus its sign determines whether the 
motion is decaying or growing; the real part represents the 
damped frequency of oscillation. 

The equations of motion for the typical section (eqs. (9a) 
and (9b)) with structural damping can be written in matrix 
form for the s th blade as 


where x a = SJmb is the distance between the elastic axis and 
center of mass in semichord units; r a = (/ a /mb 2 ) 1/2 is the 
radius of gyration about the elastic axis in semichord units; 
^ and £ a are the damping ratios; b is the airfoil semichord; 
a s is the pitching displacement; coy, = (AVm) 1/2 is the 
uncoupled natural frequency for bending; o) a = ( K a /I a ) 1/2 is 
the uncoupled natural frequency for torsion; /* and f a are the 
aerodynamic loads and s varies between 0 and jV-1. 

It is assumed that the motion of the blades is harmonic in 
time with a frequency o>and is given by 



Note that the motion has been represented as the sum of 
contributions from each interblade phase angle mode in 
which each blade has an amplitude h ar lb, a an and the phase 
angle between adjacent blades is 


G r - 2nr / N r = 0, 1, 1 (16) 


[**M} + M*} + M*} 

- {fs} = {/£} + {fas} 


The corresponding aerodynamic forces can be written in 
terms of the complex-valued unsteady aerodynamic coeffi- 
(14) cients /**,/«*, lha> laa> ^wh > an ^ Iwa : 


or 
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where p s = m 5 hrpoob 2 is the mass ratio of the blade. Thus, 
the equations for the 5 th blade are 


-a)‘ 


1 

x as 

Ml- + i® 


2 

_ x as 

r cts _ 

l a os J 


2o) hsChs 0 
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Rearranging gives 
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(17) 


To obtain the eigenvalue problem in the standard form, the 
damping terms are approximated as 


where 
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and cq 0 = reference frequency. 

To proceed further, the equations for all the N blades on 
the disk must be considered. For the assumed harmonic blade 
motion, write 
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where 


$ 

o 


© 

-> 

«o 0 


a a 0 

h oX /b 


h ax !b 

°Ol 

{y} = . 

«ai 

h o N -\ /b 


b aN-X /b 

a oN-\ , 


. tt aN - 1 , 
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Using this relation, we obtain 
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For a stability calculation (flutter), the motion-independent 
forces are set to zero, and the eigenvalue problem is obtained 
in the standard form: 


= MW + {AD} (21) 


[FMEfi]] { Y } = {0} (23) 


where 


The solution of the eigenvalue problem (23) results in 2N 
complex eigenvalues of the form 


[A/] = 


K] 


["i] 
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i— = -?= = H + iv (24) 

a>o VI 

The real part of the eigenvalue Ji represents the damping- 
ratio, and the imaginary part represents the damped fre- 
quency; flutter occurs if JT> 0 for any of the eigenvalues. 

The blade aeroelastic response induced by wakes is calcu- 
lated from equation (22) as 


{y} = -[[/’MterW (25) 




The amplitude of each blade is obtained by substituting 
equation (25) into equation (20). 
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[A) = 
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Aeroelastic Analysis for a Tuned Cascade 

For a tuned cascade (or rotor) in which all the blades are 
identical, the foregoing analysis can be simplified consider- 
ably. In this case, the aeroelastic modes consist of individual 
blades vibrating with equal amplitudes with a fixed 
interblade phase angle between adjacent blades. Hence, for 
this problem, the motion of the typical blade is written as 


Finally, after rearranging, the equations can be written as 


[['’l-MGlKn = ~{AD) (22) 



J 


j^J e «o V <V 


(26) 
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Thus, the equation for the blade becomes 


{M s ]^ lb \e i{b},+a ' s) + 





= [\]^ rlb y (m,+arS) + {AD r }e ,(b),+a ' s> (27) 


Since the blades are identical, the same equation is 
obtained for each blade. Thus, no additional information can 
be obtained by assembling the equations for all the blades on 
the disk as was done for the general mistuned system. 
Instead, equation (27) is solved for N different values of the 
interblade phase angle given by equation (16). As before, the 
equations for the forced-response problem are obtained by 
setting the motion-dependent forces to zero; the equations for 
the flutter problem are obtained by setting the motion- 
independent forces to zero. 

For the stability calculation, the equation can be simplified 
as 


([/> r ]-A[/]){y} = {0} (28) 


where 


fr] = 


+ jhhr 
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(j har 
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Stability Calculation 

The aerodynamic coefficients are calculated before the 
eigenvalue problem can be set up and solved. Since the 
unsteady aerodynamic coefficients depend on the frequency 
of oscillation, it is necessary to assume a frequency ex) 
(reduced frequency of blade vibration k c ) in advance to be 
able to calculate the aerodynamic coefficients. In actual cal- 
culations, the aerodynamic coefficients are functions of the 
inlet Mach number A/**, and the interblade phase angle a r , in 
addition to cascade geometric parameters. Either of the fol- 
lowing two procedures can be used in flutter calculations. In 
the first procedure, a value of the inlet Mach number is 
assumed, and the reduced frequency is varied until the real 
part of one of the eigenvalues Ji becomes zero while the real 
parts of the remaining eigenvalues are negative. The assumed 
flutter-reduced frequency k c f and the calculated flutter fre- 
quency vj- are both based on (Of. Thus, these two can be com- 
bined to eliminate a)/, and the flutter speed is obtained, 
namely, Vf = v J c a)jik c f. Since the inlet Mach number is 
assumed to be known, this flutter speed gives the inlet condi- 
tion (speed of sound a*,) at which the rotor will be neutrally 
stable at the assumed Mach number. The first procedure can 
be repeated to obtain a plot of flutter speed versus Mach 
number. Knowing the operating conditions, it is possible to 
determine whether flutter will occur within the operating 
region and if so, the Mach number and frequency at flutter. In 
the second procedure, the inlet conditions (Afo© and tf©©) are 
assumed to be known. A value of reduced frequency is 
assumed and the eigenvalue problem is solved. The values of 
frequency « calculated from the reduced frequency k c and 
from the imaginary part of the eigenvalue v are compared. If 
they do not agree, a new value of reduced frequency is 
assumed and the calculations repeated until the frequencies 
match. The sign of the real part of the eigenvalue fu is then 
used to decide if flutter will occur. It should be noted that the 
analysis of a mistuned cascade requires the solution of the 
equations for all phase angles at one time for a given reduced 
frequency. This is very time consuming; a procedure is being 
developed to reduce this time (ref. 66). 


where the subscript 5 identifying the blade has been dropped 
and the reference frequency cu 0 has been chosen to be equal 
to the torsional frequency <%. 

The solution of the foregoing eigenvalue problem results 
in two complex eigenvalues of the form Ji + iv (as discussed 
in Aeroelastic Analysis for Mistuned Cascades), and flutter 
occurs if Jt> 0. For the tuned rotor, the stability of each 
phase angle mode is examined separately. Hence, the 
interblade phase angle is fixed at one of the values given 
by equation (16), and the 2X2 eigenvalue problem is solved. 
The value of the interblade phase angle is then changed, and 
the procedure is repeated for each of the N permissible 
values. The critical phase angle is identified as the one which 
results in the lowest flutter speed. 
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